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CN 'Abstract 

(^ |We present a new model of scattering the plane TE-polarized light wave on an uniform dielectric layer. This wave is shown 

CN to split uniquely into two causally evolved components to describe alternative subprocesses (transmission and reflection) in all 

;_( spatial regions. Either component has one incoming and one outgoing waves, joined at the midpoints of the layer with keeping the 

M-tcontinuity of the complex-valued electrical field and the energy current density. This model, unlike the conventional one, predicts 

"^ |a subluminal energy transfer through the layer in the regime of a frustrated total internal refection (FTIR). 

^~~^ -Keywords: Hartman paradox, tunneling velocity, energy, layered structure, FTIR 



^ "1. Introduction 

•1-H . 

•*-^ , Scattering the plane monochromatic light wave on an uni- 

hr^form dielectric layer is one of the simplest scattering problems 

• °in classical electrodynamics, and, at first glance, its well known 

fj model presented in 1 1] (further named the conventional electro- 

• '-^ dynamical model of tunneling (CEMT)) represents its maxi- 
j>^mally complete and consistent description. However, studying 
f~| "the temporal aspects of this process shows that this is not the 
Q^case. Like the conventional quantum-mechanical model of tun- 

' ^nehng (CQMT) USHBIll to predict the anomalously short 

p\j 'tunneling time for a particle tunneling through an opaque poten- 



> 

in 



o 



"tial barrier, the CEMT does it for the light wave to propagate, 
■under the condition of frustrated total internal reflection (FTIR), 
'through the air gap in the structure with doubled prisms |l7|] - 
'the Hartman effect. 

In quantum mechanics the origin of this paradoxical eff'ect 
■and all difficulties to arise in studying the temporal aspects of 
'tunneling is usually associated (see |6]) with the absence in 
T— I -quantum theory of the Hermitian time operator However, in 
L| "classical electrodynamics this explanation is clearly irrelevant 
• ^H and, thus, it does not reflect the real origin of this effect - this 
^% origin must be common for both these theories. 
jd ■ Some authors (see, e.g., |8, 9]) consider that the tunneling 
■ - - "velocity concept applies only to localized wave packets, and 
the above paradoxical effect is but the result of ignoring this 
restriction. Others (see. e.g., |9]), following Sommerfeld, nar- 
row the domain of applicability of this concept to discontinuous 
wave packets. They treat the tunneling velocity as the signal (or 
information) velocity to describe the propagation of the wave- 
packet's discontinuities. 

However, the last two explanations are unacceptable too. 
Firstly, in the tunneling phenomenon ". . . an incoming peak or 
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centroid does not, in any obvious physically causative sense, 
turns into an outgoing peak or centroid. . . " ||4|, |5|]. And this 
concerns not only the wave-packet's centroid but also its dis- 
continuities, if they exist. Secondly, physics deals with the en- 
ergy or matter transfer, rather than with the information trans- 
fer Thus, above all a physicist has to take into account that the 
plane light wave transfers energy, whether it propagates in the 
infinite uniform transparent medium or tunnels through a lay- 
ered dissipative structure. Therefore, the basic velocity concept 
to characterize the energy transfer - the velocity of energy flow 
- must be applicable to tunneling without restrictions. 

Of course, this needs the individual description of the trans- 
mitted and reverberative components of the original light wave 
in all spatial regions. However, just this is beyond the scope of 
the CEMT. This model does not allow any division of the inci- 
dent light wave into these to-be-transmitted and to-be-reflected 
components, what makes impossible a correct definition of the 
velocity of transferring the electromagnetic energy through the 
layer, within the CEMT. And just this feature of this model is 
the origin of all paradoxes to surround tunneling at present. 

Thus, instead of producing new and new concepts of the tun- 
neling velocity and tunneling time, it is necessary to revise the 
CEMT in order to fulfill this requirement, and to apply the fun- 
damental concept of the energy velocity to studying the tem- 
poral aspects of tunneling. For a quantum tunneling this pro- 
gramm was realized in the approach lldllin (see also il2Lll3ll '). 
The aim of this paper is to adapt it to scattering the plane TE- 
polarized light wave on an uniform dielectric layer 

2. Setting the problem 

Let us consider two uniform nonmagnetic (//=!) media of 
dielectric permittivities eo and e: the medium of refractive index 
n in - Vf) fills the interval [0, d} on the axis OZ, and the back- 
ground medium of refractive index no ("o - V^ fiH^ the spa- 
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tial regions laying outside this interval; n,n()> \\n + hq. Both 
these media are assumed to be transparent and non-dispersive. 

Let the plane light TE wave fall from the left on the interface 
z = 0, provided that its wave vector lays in the plane YZ and the 
angle between this vector and the axis OZ is 6. In this case only 
one projection, E^, of an electrical field and two projections, 
Hy and H,, of a magnetic field are nonzero. To exploit the anal- 
ogy between the optical and quantum tunneling problems, it is 
suitable to write down these quantities to obey the same wave 
equation in the complex form (see UJ). 

Since the structure investigated is nonuniform only in the z- 
direction, we have 

£, = U(z)e'^, Hy = V(z)e'^, H, = W(z)e'^; 

X — knoyj — (jjt, no,y — riQ sin 6, k — co/c; 

c is the speed of light in vacuum. When these complex solutions 
are known the above searched-for (real) projections of electrical 
and magnetic fields are simply %{Ue'^), %{Ve''^) and 'K{We'^). 
For nonlinear characteristics - the energy density w and Poynt- 
ing vector S - of the TE wave to propagate in the medium of 
dielectric permittivity e we have 



we have also the conservation law 



w^'Uw^'\ 



S = S*"* + S 



iW 



w 



,W 



Ion L\ / J 



(1) 



(2) 



5f = 5» = 0; 



5(0) = -J-5l(f/*W); 5<_°' 



— %{U*V); (3) 

on 



Stt ^ > ' 8n ^ > 

Since the functions V and W are connected to U by the rela- 
tions (see yj]) 



V{z) = -iU'{z)/k, W(z) = -U{z)noj 



(4) 



solving the problem is reduced to finding the function Uiz); 
hereinafter the prime denotes the derivative on z- In partic- 
ular, outside and inside the interval [0,d], the initial (three- 
dimensional) wave equation for Ej, is reduced, respectively, to 
the one-dimensional equations for the function U(z), 



U" + k^nf^M = 0, U" + k\n^ - nl^)U = 0; 



(5) 



where no,z - "o cos 0. At the interfaces z = and z - d the 
function Uiz) and its first derivative U'(z) must be continuous. 
This follows from the boundary conditions for the tangential 
projections E^, Hy and for orthogonal projection H~, as well as 
from the relations (|4]i. 

Note that in the case of a plane (monochromatic) light wave, 
apart from the conservation law for the energy of electromag- 
netic field, which follows from the continuity equation 

dw 



dt 



-H VS = 0, 



5(") = VtiWU')^ const. 



ink 



(6) 



.(0) 



5 J, is the analog to the probability current density in the one- 
dimensional quantum stationary scattering problem. 

Let us write down the solutions to Eqs. (|5]), making use of 
the notations of the paper il2ll . In the region z < 0, there are 
incident and reflected waves 

Uiz) = exp(/feno,zZ) + boutik) expi-ikno^^z); (7) 

in the region z > d there is a transmitted wave 

Uiz) = ao,„ik) exp[ikno^,iz - d)] ; (8) 

inside the layer, for Q < z < d, 

Uiz) = AfG, ix-zc;k) + BfGiiz - Zc; k); (9) 



«-'-2^g* p*j' ^""•- 2\q* P*^' 



Af — Cloutj Of — Oout', Zc — ~', 

K K 2 



(10) 



Q = [G\iz) + ikG,iz)l^,^. ; P = [G'2iz) + '^G2(z)],=,^. ; 
if no,y ^ n, then 



G\ - sin(A-z), G2 - cosiKz), k - kJrf- - n^ ; 
in the case of FTIR, i.e., when no,. > n 

Gi - sinh(^z), G2 - cosh(fe), li - k Jn^ ^, - rP-. 

Here \aout?' - T is the transmission coeflicient, \bout\^ - Ris the 
reflection coefficient; T + R - I. 

Taking into account the relations (|2|i and ^ for w*"\ Sy and 
si ,we obtain 



w'">iz) = 



16n 



Sf^ = 



{n'-nl)\Uiz)fJ^ 

\U\\ 5(°U-^3(t/*t/')=^: 

Stt ~ 8nk 8n 



criQ, 



(11) 



3. The TE wave as the superposition of two components, 
transmitted and reverberative 

As in the quantum case (see flldlllh ). for any value of k there 
is a unique pair of functions Utriz) and Urejiz) which obey the 
equation 

Utriz) + Urefiz) = Uiz) 

as well as possess the following properties: (a) either function 
unlike Uiz) has one outgoing and one incoming wave; (b) the 
outgoing wave of Utriz) coincides with the transmitted wave, 
and that of Urefiz) coincides with the reflected one; (c) the 
incoming wave of either wave function is causally connected 
at the plane z = Zc to the corresponding outgoing one - the 



complex-valued functions Utr(z) and Uref(z) as well as the cor- 
responding energy current densities are continuous at this plane 
(but the first derivative of either function is discontinuous here). 
For z < 



for < z < Zc, 

U,r(z) = D,rGiiz - Zc\ k) + BfGiiz - Zc\ k) 

Uref(z) = DrefGiiz - Zc', k); 

for z> Zc 

U,r(z)=U{z), Uref = 0; 



t\(^out "out)' 



ref 



bguri'^out + bout)- 



(12) 
(13) 

(14) 

(15) 
(16) 



As is seen, the found TE waves to describe the transmitted 
and reflected TE components in all spatial regions possess the 
following properties. Firstly, not only A" + A™ , = 1, but also 
|Af"|^ -I- |AJ" ,p - 1. Secondly, the reflected TE component does 
not cross the plane z = Zc- Thirdly, despite the fact that the 
derivative U',Xz) is discontinuous at the plane z = Zc, its absolute 
value \Uii.{z)\ is continuous here, because 



\U,r(Zc-z)\ = \U,r(z-Zc)\, 
\K(Z,-Z)\^\K(Z-Zc)\. 



(17) 



These relations follow from the equality %(DtrB*,) = %{AfBy) 
which follows, in its turn, from Exps. ( fTOl l and ( fl5b . 

All this means that for the transmitted TE component the real 
fields E'^ and H[!', as well as the energy density w"' and the 
Poynting vector S"^, are continuous at the plane z = Zc'- 



s"' = (o,s';,s'!y, s'; = ^|f/„f 



"""^^'-Xe-. 



Stt 

[r?+nl)i\U,r{zf 



S'! = 5f ; 

if/;(z)r 



(18) 



k^ 



But the real projection //"' is discontinuous at this plane. 

So, the light wave spUts uniquely into two causally evolving 
components, transmitted and reverberative. Now, when each 
component is known in all spatial regions, we can proceed to 
studying the temporal aspects of this scattering process. 

4. Transmission and reflection dweU times 

The velocity \"'{z) (v"' - (0, v'f , v'/)) of the light component 
to propagate through the layer is introduced, in our approach, 
as the ratio of the Poynting vector to the energy density at the 
point z (see Exps. ^y. v"'(z) = S"'(z)/w"'(z). Thus, the time 
T^ to describe the duration of the transmission of the light wave 
through the layer can be defined as follows 



'°=r^=s?r»'"«* 



(19) 



hereinafter, this quantity will be referred to as the transmission 
dwell time. 

Note that flie CEMT deals wifli the Buttiker dwell time tq. 



1 r (0) "-"0 ' 

TD = -;— W^ '(Z)dz, line = -^, 

line Jo OJT 



which does not distinguish between the to-be-transmitted and 
to-be-reflected components of the incident light wave. This 
fully concerns also the so called "group tunneling time" tgs (see 
il4ll and references therein). 



1 



w^'^\z)dz. 



(20) 



Despite different normalization, this quantity has no relation to 
tunneling, because, like t/j, it is defined via w^^\z). 

Then, with taking into account (I?])-® and (fT2l i-(fT6]l. for the 
transmission time we obtain the foUowing expressions. For 
no,v < n 



Tn = 



- Un^ - nl)nl ^ sin(2Kd) 



+2n (n + Hq , — Hq ^)Kd\; 
for «o,v > n (the FTIR case) 



(21) 



-[(«„- n^)r?Q sinh(2A-t/) 



Alc'crid,, ' 

— 2n (n + «() ^^ — Hq y)Kd\. 



(22) 



Figs.[T]and|2]show the dependence T'jjJTfree on kd for the case 
when one medium is vacuum, and another is glass; T/^e = die. 
As it follows from Exp. (l22l) (see also the curve 4 on Fig. |2l 
obtained for no = 1,5 and n - 1), in the case of FTIR {6 > 
41,8°) this quantity exponentially increases when kd ^> oo. 
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Figure 1: The dependence of r'^JTfree onkdfovrto = 1 andn = 1,5: (1)6 = 0°; 
(2) e = 15°; (3) d = 30°; (4) 8 = 45°. 
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Figure 2: The dependence oir'^/Tfi-ee onkdfoi no = 1,5 and« = 1; the values Figure 3: The dependence i'"/c on z for no = I and n = 1,5; kd = 10; the 

of the angle 9 are the same as for Fig.[T] values of the angle 9 are the same as for Fig.[T] 



When «oj < « the transmission velocity does not exceed c 
too (see Fig.[TJ. As it follows from (ISTT i. in the limit kd —^ co 



-^fn 



2Kkno,z 



yn^-^ 



> 1. 
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Note that t^ does not yield a full information about the veloc- 
ity of the energy propagation inside the layer, since r^ depends 
only on the z-projection v'S. Therefore a more detailed analysis 
of this question is done in Section |5] 

The reflection dwell time t^^ is defined similarly - 



ref 



rf^i 



reft \J cref '^"O.z „ 

w H,z)dz; S,' ^ "stT 



1 



lOTT 



( 1 2\ 2 KaM)\^' 



k^ 



For «() V < n, with considering Exps. ([TZb-lfTSIl, we obtain 
2no,2 "^'^'^ ~ "o,y sin(/crf) 



ref 



for noy > n 



CK («2 _ „2) COS(Kd) + n^ + Hq , - Hg 



(23) 



ref 



2ni 



Q,z 



iIq , sinh(/f(i) - n iid 



CK («2 _ „2-) cosh(/frf) - (n^ + ^2 - „2 ■) 



From Exps. (l23T l it follows that, under the conditions of FTIR, 
the function t'^ (d) saturates in the limit d — > oo. However, this 
fact does not say that we deal with the Hartman effect, because 
the reflection time depends not only on the velocity of the rever- 
berative TE-component, but also on the depth of its penetration 
into the layer. The above fact means simply that in the case of 
FTIR this depth tends to some fixed value when t/ — > oo. 



5. On the velocity of the energy transfer 

A detailed analysis of the velocity of the energy transfer for 
the light component to pass through the layer has been carried 
out by the example of the above two media - vacuum and glass. 
Figs. [3]|5]present numerical results obtained for the case when 
the interval [0, d] is filled with glass. Fig. [3] shows the func- 
tion v"'(z) = |v"(z)| within the layer. Fig. |4] displays the z- 
dependence of the angle (0 = aictan{S [^ / S l^)) to characterize 
the direction of propagation of the transmitted TE component. 




Figure 4: The dependence of the angle on z for the same parameters as for 
Fig.S 



As is seen from Figs. [3] and 31 both the functions - v"'(z) and 
0(z) - reach their maximal values on the set of points, which in- 
cludes the boundary points z = and z - d. When increases. 
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Figure 5: The functions v"'(z)lc (firm line) and v^""(z)/c (dashed line) for 6 = 
89, 1°; the values of remaining parameters are the same as for Fig. [3] 



Figure 6: The dependence of v" on z for no = 1, 5 and n 
values of the angle 6 are the same as for Fig.[T] 



1; kd = 10; the 



these functions vary more rapidly. In this case, 0(z) < 9 for 
«o < «• 

It should be noted that the refraction angle, 0oo(0), to charac- 
terize scattering the plane light wave on the semi-infinite dielec- 
tric medium has no direct relation to the process under study, 
where this dielectric fills the layer of a. finite width. For a given 
6 the function 0(z) oscillates, inside the layer, around the value 
which is approximately equal to ®oo{9)- 

In the limit — > 90", the function v"'{z) tends to zero at the 
points of minimum. This is seen from the numerical results 
presented on Fig. [3] as well as on Fig. |5] where, in addition 
to v"'(z), we show also the function v^""(zy, v^""(z) - \\^""(z)\; 
\f""(z) = S(°Vw<">(z). At the point z = this function like v"'(z) 
is discontinuous, but its discontinuity is so small in this case 
that it is unapparent on the figure. 

Fig. |5] shows explicitly the principal difference between the 
behaviour of the functions v"'(z) and v^""(z) near the interfaces 
z = and z = d. Due to ( fTTT i \"'(zc - z) - v"'(z - Zc)- However, 
the "velocity" v^""(z) to underlie the concept of the "group tun- 
neling time" ( l20t . does not obey this requirement! This fact 
presents one more argument in favor of our approach. 

Indeed, in the symmetrical structure to consist of transparent 
uniform media, all reflection symmetric points are physically 
equivalent. Thus, the tunneling velocity must represent the even 
function of z - Zc- The reverberative light component to exist 
only in the left half of the interval [0, d] must not affect the 
tunneling velocity. 

Of importance is to stress that v"'(z) = c/no outside the inter- 
val [Q,d]. Inside this interval the function v"'(z) varies. How- 
ever its values do not exceed here the limiting velocity c. For 
n > no the velocity v" (z) takes its maximal value at those points 
z where sin(/cz) - 0. This set is always nonempty, as it contains 
the boundary points z = and z = li. At any point of this set 
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Figs. |6][8]present numerical results for the case when, in the 
considered structure, glass and vacuum traded places. Figs. |6] 
and Q show, respectively, the functions v"'(z) and 0(z). Fig. [8] 
shows the functions v''(z) and vi" (z) for their comparison in 
the case n < «o_v 

Note that for no > n > no,v the velocity v"' takes maximal 
value at those points z where cos(/cz) = 0: 



v"- = v®, . = - . 



2k^noM Jk^ + k^n^ n^ 



n k^n^M^ + ^K^ + k^nl^nl^) « 



c 

^ n 



If n = no^y, then the function v" (z) reaches the maximal value 
c/n at the point z - Zc- This is the only case when the velocity 
of light in the medium, fitting the finite layer [0, d], approaches 
its velocity in the same medium, but fitting the infinite space. 

In the case of FTIR, the maximal value v"'(zc) diminishes 
when the angle 6 increases; if exceeds some critical value, the 
function v"'(z) reaches its maximal value v;„a^ at the boundary 
points z = and z - d. As regards the point Zc, in the limit 
-^ 90°, v"'{zc) - c/no (see Fig.[8]l. When no > n the inequality 
0(z) >0 holds (see Fig. Illl. 

6. Conclusion 

A new model of scattering the plane TE-polarized light wave 
on the uniform dielectric layer has been developed. It is shown 
that this wave can be uniquely presented, in all spatial regions, 
as the superposition of two causally evolved transmitted and 
reverberative components. Unlike the original light wave either 
component possesses one incoming and one outgoing waves, 
joined on the plane z = Zc with keeping the continuity of such 
quantities as the energy density, the density of the energy flow, 
the (real) electrical field and the y-th projection of the (real) 
magnetic field. The z-th projection of the (real) magnetic field 




Figure 7: The dependence of the angle on z for the same parameters as for 

Fig.ia 



for either subprocess is discontinuous at this plane (however, 
their sum is continuous here). 

Within this approach, the concept of the transmission veloc- 
ity can be consistently introduced as the velocity of the energy 
transfer, via the Poynting vector and energy density to describe 
the transmitted component. Unlike the existing CEMT's con- 
cepts, ours obeys two important requirements: it does not lead 
to the Hartman effect in the case of FTIR; it represents the even 
function of z-Zc for the structure to possess the mirror symme- 
try, with the symmetry plane z = Zc- 

Of importance is that the present model is applicable to scat- 
tering the plane electromagnetic wave on a layered structure 
with dispersion and dissipation. 
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Figure 8: The functions v"'(z)/c (firm line) and v^""(z)/c (dashed line) for I 
89, 1°; the values of remaining parameters are the same as for Fig.|6] 
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